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Abstract. Multipolar solutions of Maxwell's equations are used in many practical 
applications and are essential for the understanding of light-matter interactions at the 
fundamental level. Unlike the set of plane wave solutions of electromagnetic fields, 
the multipolar solutions do not share a standard derivation or notation. As a result, 
expressions originating from different derivations can be difficult to compare. Some of 
the derivations of the multipolar solutions do not explicitly show their relation to the 
angular momentum operators, thus hiding important properties of these solutions. In 
this article, the relation between two of the most common derivations of this set of 
solutions is explicitly shown and their relation to the angular momentum operators is 
exposed. 
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1. Introduction 

The solution to Maxwell's equations in homogeneous isotropic media is one of the most 
studied problems in undergraduate physics. Often the problem is solved in Cartesian 
coordinates and gives rise to the typical plane wave solutions of electromagnetic 
fields. Although other solutions in different coordinate systems are not so widespread, 
the multipolar fields arising from the use of spherical coordinates stand out due 
to their use in several important problems in electromagnetism [TJ. These sets of 
electromagnetic fields play a very important role in the interaction of light with matter. 
In particular, they are widely used to understand the atomic and molecular spectra 
[2j [3]. More recently, the multipolar fields have found interesting applications in the 
field of nanophotonics [4]. Probably the most successful application of this set of 
solutions is the Mie scattering problem, which studies the interaction of electromagnetic 
plane waves with spherical particles (5]. Mie theory provides one of the few analytical 
solutions of the inhomogeneous Maxwell's equations. This theory generalizes Rayleigh 
scattering, which accounts for the colour of the sky, to particles of any size. For all 
these reasons, the multipolar fields appear in almost all undergraduate textbooks and 
in a long range of specialized books. Unfortunately, the multipolar fields do not share 
the structural simplicity of the plane wave solution. This fact translates into a myriad 
of different representations for the multipolar solutions. Different books use notations 
and derivations which are so dissimilar that the study of these solutions can be difficult. 
More importantly, some of the standard derivations hide the beautiful symmetry of the 
multipolar fields and can prevent a deeper understanding of the underlying physics, i.e. 
the angular momentum of light. 

In this article we would like to present the most common derivations of the 
multipolar fields and how the different notations are related. This exercise usually 
requires a significant amount of algebra and is rather cumbersome [6]. Here, we will 
show that an insightful understanding of the role of the angular momentum operator 
in electromagnetism and the proper use of some of its properties significantly simplifies 
the calculations. With this new understanding, the multipolar solutions become much 
clearer and can be used in physical problems with a new point of view. 

The remainder of the article is organized as follows: In section [2] the problem is 
outlined, along with two common derivations of the multipolar solutions. By making 
use of the orbital angular momentum operator, the correspondence between these two 
sets of solutions is discussed in section |3j Finally, section H] concludes the work with a 
summary and outlook. 

2. Multipolar fields 

Let us start by reviewing the multipolar solutions to Maxwell's equations for source- 
free homogeneous and isotropic media. We consider monochromatic electric, S, and 
magnetic, H, fields, i.e. £ = 9fte{Eexp(— iut)} and similarly for the magnetic field 
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TL = 5ie{Hexp(— iut)}, where the vectors E and H follow Maxwell's equations and 
therefore the resulting vector wave equations 

V 2 E + fc 2 E = 0, V 2 H + A; 2 H = 0, (1) 
with k = ^/ejluj, have zero divergence 

V-E = 0, V-H = 0, (2) 
and are interrelated: 

V x E = iufjtH, V x H = -iueE. (3) 

Here e = eoe r , where eo and e r are the vacuum and relative permittivity, respectively, 
and similarly /i = /io/i r , where fiQ and fi r are the vacuum and relative permeability. 
Multipolar fields are solutions of the above equations in spherical coordinates. They 
may be chosen to be simultaneously eigenfunctions of the squared angular momentum 
of the radiation field, a component of the angular momentum on one of the coordinate 
axes (here, the z-axis is chosen), and of the parity operator. The multipolar fields are 
often used as basis functions, into which electromagnetic fields can be decomposed [7J. 
Specification of the weights in such a multipolar expansion allows the description of 
electromagnetic fields. 



2.1. Derivation of the multipolar solutions of Maxwell's equations 

One can find several expressions for the multipolar solutions of Maxwell's equations in 
the literature, but they are typically derived in two different ways. The first one uses 
a well-known recipe for deriving the solution of ([I]), ([2]) in a medium free of charges [8]. 
This recipe is applicable to different coordinate systems that meet certain conditions 
[8], and consists in using the following fields: 

M = V x (ur), N = V * M ; (4) 

k 

where u is a unit vector in the symmetric direction and r is a scalar function expressed 
in the appropiate system of coordinates. This scalar function r can be written as 
t = w(x)ip, where \ is the coordinate associated with the available symmetry and 
w(x) — depending on the problem; if) is a solution of the scalar Helmholtz 

equation in the chosen system of coordinates. If all these conditions are fulfilled, the 
vector fields M and N are solutions of Maxwell's equations showing the appropriate 
symmetries. 

In the case of the multipolar solutions, the system of coordinates is spherical, the 
unit vector points in the radial direction, u = f, and w(r) = r. This approach is 
typically found in Mie theory related literature [7J as well as in some electromagnetism 
textbooks jU |9j [ID]. We will refer to this method as the "general method," since it may 
also be used to obtain solutions in other geometries. 

Another way of deriving the multipolar solutions of Maxwell's equations uses the 
fact that these solutions have to be eigenvectors of the angular momentum operator. 
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Hence, the solution can be built from the rules of addition of the angular momentum 
operator. This construction is detailed in the books by Rose [TH [12]. In this method, 
the multipole solutions are explicitly shown to be eigenvectors of the angular momentum 
operator, a fact that is typically overlooked in many books about Mie theory. We will 
refer to this method as the "angular momentum method," as it explicitly makes use of 
the angular momentum operator properties. 

The diametrically different construction methods of the multipole solutions may 
lead to some confusion due to several factors. On the one hand, as we will see 
immediately, both methods naturally give rise to different notations, which are rather 
cumbersome to reconcile [6]. Then, one may be tempted to think that the different 
multipole notations belong to different classes of solutions when in reality they represent 
the same object. On the other hand, the derivation of the multipolar solutions without 
making use of their role as eigenvectors of the angular momentum operators hides 
some important properties of this set of solutions. For example, it is often difficult 
to understand the relation between the different resonances in Mie scattering and the 
angular momentum of the field and/or the particle. In this paper, we will try to solve 
this problem once and for all and show how the two construction methods are related. 



2.2. General method 

As mentioned before, the general method is based on (jlj). The scalar function r is one 
of the solutions of the scalar Helmholtz equation in spherical coordinates multiplied by 
the radial coordinate (r = rip), and the unit vector is oriented along the radial direction 
(u = f). Our notation closely follows that of Bohren and Huffman [7], which is very 
common in Mie theory books. Note, however, that we distinguish different modes with 
the indices M and L, instead of the m and n that are used in the reference. The choice 
for this trivial change will become clear in the following sections. Using (jlj), we define 
a set of solutions {NL eML ,~M oML ,N eML ,~N oML }, which are determined by the scalar 
functions tp e ML = cos(M(j))P£f (cos(8))z L (kr) and ip oM L = sm(M(j))P£f(cos(8))z L (kr). 
Here, as in the rest of the article, (r, 8, 0) are the spherical coordinates (0 being the 
azimuthal angle). The radial function zi{kr) can be either a spherical Bessel function 
or spherical Hankel function. Within the general method the indices M, L are non- 
negative integers and L > M. As required by the recipe of (jlj), the functions ipoML 
and ipeML fulfill the scalar Helmholtz equation in spherical coordinates. The definitions 
required to construct the complete form of the solutions {M eML ,'M oML ,'N eML ,N oML } 
can be found in the appendix. As an example, 

M oML = —j- cos(M0)P L M (cos(#))z L (^)e e - sin(M0) — { ° S{ ^ z L {kr)e^ (5) 

where e# and are the polar and azimuthal unit vectors. 

Note that there are two different sets of M and N solutions labeled with either 
o (odd) or e (even), and they differ in the azimuthal behaviour of the scalar function 
ip, oscillating as a sine or a cosine, respectively. The main difference between these 
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multipolar fields and those given by Jackson [TJ is that instead of sinusoidal functions, 
the author uses complex exponential functions. 

2.3. Angular momentum method 

As mentioned earlier, the general method does not make use of the fact that the 
multipolar solutions should have a well defined relation with the angular momentum 
operator!]] However, one can derive the multipolar solutions by directly imposing that 
they should be eigenvectors of the angular momentum operator. We will follow the 
notation by Rose [HJ [12] ■ First we start by noting that the angular momentum operator 
(J) for a vectorial field is composed of two terms: the orbital angular momentum 
L = —i(r x V) and the spin angular momentum S, i.e. J = L + S. In this article, we use 
dimensionless angular momentum operators, as in the references cited jl] [TTJ [12] . The 
angular momentum operator so defined is the generator of rotations for the vectorial 
electromagnetic field. It is also related to the density of electromagnetic angular 
momentum of the field (i.e. m = Jie{r x E x H*}). The interested reader can find 
more details about these subtleties in the given literature. 

We should now impose that the multipolar vector fields are eigenvectors of J 2 and 
J z . We can start constructing these solutions by composing the eigenfunctions of both 
L and S. The eigenfunctions of L 2 and L z are the spherical harmonic scalar functions 

(Y L M (e, <!>))■. 

\?Yl l = L(L + 1)F L M , L z Yl 4 = MY™. (6) 

Note that the above equations can be considered as the definition for the spherical 
harmonics as they are a set of partial differential equations. On the other hand, 
the eigenfunctions of S 2 and S z are the vectors £ +1 = — (x + zy)/v2, £ = z, and 
= iy)/V2: 

S%. = 2^, S z ^ = (7) 

We then construct the multipolar vector functions Tjlm by combining the scalar 
eigenfunctions (Y z * /+/1 ) and the vector eigenfunctions (£_„), using the rules of angular 
momentum addition. A simple product of the scalar and vector eigenfunctions, 
Y^ +ti £_ fl , remains an eigenfunction of the orbital and spin angular momentum 
operators. However, it does not have the required properties with regard to the squared 
total angular momentum operator, J 2 . In order to create eigenfunctions of the squared 
total angular momentum operator, one has to combine the separate eigenfunctions with 
the aid of the Clebsch-Gordan coefficients: 

At=+1 

Tj LM = C ( 1LJ ~> M + ») Y L^{0, 0)£-m> (8) 

| Of course, the angular momentum content of the multipolar solutions can always be calculated and 
gives the expected result p]. 
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where C(abc;d,e) denotes the well-known Clebsch-Gordan coefficients, which we have 
explicitly written in the appendix. Note that T JLM is an eigenfunction of J z , because the 
three terms in the sum above (/i = —1, 0, +1) combine eigenfunctions with L z = M + fi 
and S z = — always summing J z = L z + S z = M. Hence, in this notation M has 
become the eigenvalue of the z-component of the total angular momentum operator. 
Furthermore, the weights of the three components in the sum are the Clebsch-Gordan 
coefficients, which are exactly calculated so that the total angular momentum fulfills 
3 2 Tjlm = J (J + 1)Tjxm- In summary, the vector functions Tjlm are eigenvectors of 
both J 2 and J z , but also of L 2 and S 2 . They are often referred to as vector spherical 
harmonics and are also used by Jackson p], even though the derivation in this reference 
is different. A further noteworthy property of these vector spherical harmonics is that 
they are orthonormal on the surface of a unit sphere. 

The electric and magnetic multipolar solutions can be derived from the vector 
spherical harmonics in (JSJ). The addition rules for angular momenta impose that for a 
fixed J, the orbital angular momentum L can only take three values: J— I, J, and J+l. 
To proceed, we can make use of the fact that the multipolar fields are also required to be 
eigenvectors of the parity operator, P. This means that they are invariant under point 
inversions (r — > — r), except for possibly an overall phase difference. The parity of Tjlm 
is determined by that of the scalar spherical harmonics Y^ +fl , which is (— 1) L . Then, 
to ensure that the fields are eigenfunctions of the parity operator, only terms with equal 
parity can be combined, i.e. the solution with L = J must remain separate, while those 
with L = J ± 1 can be mixed. The multipoles are classified as electric or magnetic poles 
according to their parity: If A has parity (— 1) J , then it is a 2 J pole magnetic field, and 
in the case of A with parity (— 1) J+1 it is a 2 J pole electric field. For this reason, the 
magnetic multipole vector potential can only be composed from a term with J = L. 
The vector functions used are therefore 

TW = Yl C ( 1LL > M + a0^ m+m (^ (9) 

Then, in order to obtain the vector potential for the magnetic multipole field, Tllm 
has to be multiplied by a spherical Bessel or Hankel function: 

A LM (m) = C L {m)z L (kr)T LLM . (10) 

Cx,(m) is a normalization factor, which is set to -1 in [12j. Now it is very easy to obtain 
the electric multipoles by using (jHJ) or the symmetries between the electric and magnetic 
fields. 

By construction, the electromagnetic fields with this shape are eigenvectors of the 
angular momentum operators. Now, having followed the two derivations, it is unclear 
how the fields obtained by the general method are related to the ones from the angular 
momentum method. The final expressions for the fields from the general method (see 
for example ([5])), which are typically given in terms of spherical unit vectors, have a 
dissimilar presentation to those from the angular momentum method (e.g. (TlOl) ). As an 
illustration of the difference in appearance between the two sets of solutions, the squared 
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modulus of each component of the field {L, M} = {3, 3} on a spherical surface is shown 
in figured] Comparing the corresponding field components from the two formalisms, it is 
evident that there is a difference in both azimuthal dependence and magnitude between 
the fields. In order to determine the relationship between the two sets of solutions, one 
could use a brute force approach to calculate the angular momentum of the fields of 
the general method and to calculate the conversion factors from one field to the other 
[6]. Here, we will present a different method by using the properties of the angular 
momentum operators. 




Figure 1. (Colour online) Squared modulus of the x- ((a),(d)), y- ((b), (e)), and z- 
((c), (f)) component of the multipole field with indices {3, 3} on a spherical surface in 
each formalism. Top row: M D 33 from the general method. Bottom row: £33(771) from 
the angular momentum method. 



3. Correspondence between fields 

3.1. Properties of the angular momentum operator 

We will start by reviewing some of the properties of the angular momentum operator, 
which we will use later. The angular momentum operator is a vector differential operator 
that can be written in the form: 

L = -i(rxV) (11) 
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= L x ± + L y y + L z z (12) 
= -L+x^ + L Z £ -L_!£ +1 . (13) 

Note that ( Hip is very similar to the quantum mechanical angular momentum operator, 
normalized by h, i.e. L = r x p/h~. The linear momentum is written in its differential 
form p = —iHV. Equation ( !T3|) can be easily verified with the form of the vectors £^ 
and using the so-called ladder operators: 

L+i = - -j=(L x + iL y ) 

(15) 

From the theory of angular momentum it can be shown that, when applied to 
spherical harmonics, the ladder operators have the following effect: 



:(L x -iL y ). (14) 



j y m = i (L-M)(L + M + l) M+] 
\ i . y 2 / 

s/ / (L + M)(L-M+l) v , ; _ j 



L iY m = y v- ■ ~, v - ~ ) (i6) 

and indeed they can be used to obtain the relations between spherical harmonics of 
different orders and their explicit shape. 

The final property of the angular momentum operator that we will make use of is 
that as the generator of rotations, the operator does not affect the radial components 
of the fields. This property is very easy to verify by writing the angular momentum in 
radial coordinates or by observing that (fTEj) operate on the spherical harmonics, which 
are purely angular functions (they contain no radial component). 



3.2. Rewriting the general method 

Now we can turn to the general method again and observe that the fields M can be 
rewritten in the following form: 

M oML = V x (r^Wi) = -r x (Vt/Wl) = -iL^ oML , (17) 

and similarly for the even fields, i.e. M £ ml = —iL^eML- The second equality in (fTTl) 
can be obtained from the commutation relations between r and p or directly by checking 
terms of the kind yd/dz — zd/dy = —(d/dyz — d/dzy). 

Finally, from the explicit expression of the spherical harmonics, which is given in 
the appendix, it can be shown that ip e ML + iipoML oc Yff ZL(kr) and ip e ML — iipoML oc 
Y L M zi(kr). Using this and the properties of the angular momentum operator given 
above, we can combine the fields M in the following way: 



M eML ±iM oML = -i(TlY M \l 1 [ L L + ^ z L (kr)LY L ^( 
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-iz L (kr)(Tl) 



(M) 



4tt (L + M)\ 



2L+1(L-M)\ 
+ L z £, - L^ +1 )Y± M (e, 0), (18) 

which now can be cast as a sum of different spherical harmonics by using ([16]) : 



M eML ± iM oML = ~ iz L {kr)(Tl) 



M I 47T (L + M)\ 

2L + 1 (L-M)\ 



(L±M)(L^M + \) Y±M _ 



(19) 



Thus, we have been able to easily write the M fields of the general method as a 
superposition of products of angular momentum eigenstates, similarly as in the angular 
momentum method. Now it is just a matter of comparing the factors in ( fl9l) . obtained 

P and (IA.11D - 0A.13P and noting that they 
i.e. 



through the use of ladder operators, with 
are the same, except for a common factor! 



M eML ± iMoML = K^) M ff 2 ^ f L + _ M ] f L(kr)T LL±M , (20) 

which gives us the desired relationship between the two methods. Let us consider (1201) 
in more detail. The constant factor on the right hand side is due to the difference in 
magnitude between the two sets of solutions. Based on properties of T LL ± M established 
in section 12. 3[ (|20|) clearly shows that M eML ± iM. oML is also an eigenfunction of the 
operators J 2 , J z , and P. Conversely, note that the modes resulting from the general 
method, i.e. {M eML ,~M oML ,~N eML ,N oML } are eigenvectors of J 2 and P, but not of 
J z . This is reflected by the difference in azimuthal dependence between the two sets of 
solutions, visible in figure (U which is accounted for by the use a linear combination of M 
fields in ( l20l) . Note that the need for this linear combination stems from the particular 
choice of ^-dependent functions during the solution of the scalar Helmholtz equation in 
the derivation of the fields from the general method. If complex exponential functions 
were chosen instead of the sinusoidal functions, the fields derived by the two methods 
would be directly proportional to each other. 

We could similarly derive another relation for the N fields, but it is easier to use 
that N = (V x M)/fc and H = — i(V x E)//c[Jj| as will be evident from the conversion 
tables below. 

§ That they are the same is not a lucky coincidence. The Clebsch-Gordan coefficients are precisely 
derived by using the angular momentum ladder operators. 
II In [11] and [12] Gaussian units are used. 
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3.3. Conversion tables 

We can now present the whole conversion tables between the electric and magnetic 
multipoles from Rose's books (denoted by (e) and (m), respectively) and the M and N 
fields from the general method used in Mie theory. It is important to emphasize that 
L(L + 1) and M represent the eigenvalues of the square and z-component of the total 
angular momentum operator, J. The vector potential can then be written as: 



AlmM = n-77!#^ V^yS(M e iM|x + zsign(M)M 



(21) 



A -( g ) = ( ^^ ^»( N ^ + 1 ^ w*w 

And for the electric and magnetic fields: 
E IM (m) = -H LM (e) = —J^^ ^mT^^^ + , sign (M)M o|M|I ) 



3.^. Application to a plane wave 

As previously mentioned, the plane wave solution is very commonly used in 
electromagnetism. A physical problem that ties together the use of multipolar fields 
and plane waves is Mie scattering. The incident field is a plane wave, but due to the 
symmetry of the spherical scatterer it is convenient to treat the problem using multipole 
expansions, which necessitates the decomposition of the incident plane wave in terms of 
the multipolar fields. In the following, we demonstrate the use of the above conversion 
tables for the case of a plane wave. To convert between the two representations, we 
can use (12TT) on the plane wave decomposition given in [TT] to obtain the corresponding 
expansion in terms of {M. eML ,~M oML ,N eML ,~N oML }. In [TT], the plane wave expansion 
is given as 

oo 

A pw = v / 27^2 L v / (2^ + l)(A Lp (m)+^AL P (e)), (23) 

L=l 

where the vector potential for a plane wave is A pw = u p e tkz and u p = ^(x + ipy), 
p = ±1. Choosing ^(ui + U-i) to obtain a plane wave with polarization in the x 
direction, and applying ( l2Tj) yields 

oo 

A pw = ^J2 tL V( 2L + l)(A L1 (m) + iA L1 (e) + A x _i(m) - zA L _ x (e)) 

L=l 

oo 



tl V L (L + 1) V 47r ( L + l ) L 

[-i(M elL + iM olL ) - i(N elL + iN olL ) + z(M elL - iM olL ) - i(N elL 
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The difference between these expansions is illustrated in figure [2J which depicts the 
weights of the first few basis functions with index {L, ±1} in the representation from 
the angular momentum method, and {1, L] in that from the general method. 



(a) 
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Figure 2. (Colour online) Complex weights of the first basis functions in the multipolc 
expansion of a plane wave, (a) Vector potential using Rose's notation, where weights 
multiply (Aii(m) + iAi,i(e) + Al_i(to) — iAz,_i(e)). (b) Electric field in the notation 
from [7], where weights multiply (M g il — zNeii). 



4. Conclusion 

In this article we have compared different notations of the expressions for multipolar 
modes in electromagnetism. We have exposed their symmetry properties, even in the 
cases where they were hidden in the formalism. As a result, we have given compact 
explicit expressions that relate the different notations. All this has been possible by 
making use of the angular momentum operator and its properties. The relation between 
the multipolar electromagnetic modes and the angular momentum is typically overlooked 
in Mie theory, which may lead to some confusion. By presenting some of the most 
commonly used expressions for multipolar fields on equal footing, we hope to have lifted 
this possible confusion. 
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Appendix 

For the convenience of the reader, here we will write the expressions of some of the 
functions we have used above. 

Associated Legendre functions 
The associated Legendre functions used in [7J are defined as follows: 

( 1 _ ~2\Af/2 jM+L 

M 4^ - 1)L - (A - 1) 
The form used in spherical harmonics has x = cos(9), which is 

This can be seen from page 90 of [TJ, which provides the definition in terms of the 
ordinary Legendre polynomial. However, another common definition of the associated 
Legendre polynomial, e.g. see pQ, has an additional factor of (— 1) . Therefore, care 
must be taken when comparing different references. 

Spherical harmonics 

The most common definition of spherical harmonics, also used in [TT] (p. 241), is 

Using the definition of associated Legendre polynomials given by (1A.2|) . this means 

that 

Y L M (9,4>) = (-1) M ^/ 2 ^ 1 f L ~J ^Pl (coB(g))e^. (A.4) 
Based on this, 

cos(M0)P^(cos(0)) + zsin(M0)P* f (cos(£)) = e^P^(cos(0)) 

cos(M0)Pf (cos(0)) -isin(M0)Pf (cos(0)) = e~ iMtt> P^ (cos(0)) 

= e-^P L ^(cos(0))(-l)^^±^| 

/ 4tt (L + M)!^ „ M ^ 
2L + 1(L-M)! ^ (A ' 6) 
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Spherical Bessel functions and Hankel functions 
The two spherical Bessel functions relevant for multipole expansions are defined in terms 
of common Bessel functions as 



h{kr) = x /^J L+ i/2(kr) (A.7) 



y L (kr) = ^Y L+1/2 (kr). (A.8) 

Their linear combination is used to construct the spherical Hankel functions 

ftW (Ax) =j L (kr) + iy L (kr) (A. 9) 

hf (kr) = j L (kr)~ iy L {kr) . (A. 10) 



Clebsch-Gordan coefficients 



C(1Li;0 ' M) (A - 12) 

C{\LL; 1, M — 1) - f + (A,3, 
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